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Abstract We present a detailed investigation of the behavior of the nonlinear q- 
voter model for opinion dynamics. At the mean-field level we derive analytically, 
for any value of the number q of agents involved in the elementary update, the 
phase diagram, the exit probability and the consensus time at the transition point. 
The mean-field formalism is extended to the case that the interaction pattern is 
given by generic heterogeneous networks. We finally discuss the case of random 
regular networks and compare analytical results with simulations. 
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1 Introduction 

Ordering dynamics is a classical subject of nonequilibrium statistical physics, re- 
garding phenomena such as magnetic ordering or phase-separation of mixtures E] 
|2T1 . Also in the biological and social domains the question of how order emerges 
and evolves over time is of fundamental interest. This has recently led to the intro- 
duction of a wealth of models for ordering dynamics in interdisciplinary settings, 
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inspired by the physical paradigms but exhibiting a diversity of specific distinctive 
features ifTTl . In this context the voter model [14 20] plays a prominent role for 
several reasons: it is extremely simple; it has a natural interpretation in several 
fields (ecology, genetics, opinion dynamics); it can be exactly solved in regular 
lattices for any number of spatial dimensions II22II21I . 

The voter model is defined as follows: A population of agents is placed on 
a regular lattice or a graph of size N. Each agent is endowed with a binary spin 
variable, representing two opposed opinions, and taking values a = ±1. At each 
time step t, an agent i, together with one nearest neighbor j, are selected at random, 
and the state of the system is updated as a, := Gj, the first agent copying the 
opinion of its neighbor. Time is updated t — > t + 1/N. In this way, the flipping 
probability, i.e. the probability that a spin surrounded by a fraction x of spins in 
the opposite state flips, is f{x) = x. Starting from a disordered initial state, this 
dynamics leads in finite systems to a uniform state with all individuals sharing 
the same opinion (the so-called consensus). From the social dynamics point of 
view, interest is mainly focused on the exit probability E(x) and the consensus 
time Tn(x), defined as the probability that the final state corresponds to all agents 
in the state +1 and the average time needed to reach consensus in a system of 
size TV, respectively, when starting from a homogeneous initial condition with a 
fraction x of agents in state +1 Hill . On regular lattices voter dynamics conserves 
on average the global magnetization m = Y4L1 an d this implies E(x) = x. 

On the other hand, for d > 2 J3 the consensus time takes the form 

T N {x) = -N eS \xhi{x) + (\ -x)ln(l -*)] , (1) 

with N eff ~ N 2 in d = 1, N eS ~ NlogN md = 2, and N eS ~ N in d > 2 (the latter 
being also the mean-field result) 1211 . 

In the social context voter dynamics becomes all the more relevant when it is 
considered on top of complex networks, which are more realistic representations 
of social contact or relationship patterns [3 25 |. The analysis on these substrates 
reveals nontrivial differences with respect to lattices. For example, now the order 
in which interacting individuals are selected matters H3 Oil 101 . Moreover, relevant 
quantities such as the exit probability and the consensus time turn out to depend 
on the heterogeneity of the contact pattern, as measured by the degree distribution 
P(k) [ 29 28 27 ]. On complex networks the main conserved quantity is not the total 
magnetization, but a connectivity weighted magnetization, ft). Thus, it is possible 
to show that the exit probability is equal to this conserved quantity, while the 
consensus time takes the form [29] 

ZV(fl>) = -Nat [(Bln(ffl) + (1 - ft))ln(l - ft))] , (2) 

where now N eS = N(k) 2 /(k 2 ), with (k n ) = T, k k n P{k) the n-th moment of the de- 
gree distribution. 

In the case of regular lattices, the simple behavior of the voter model has been 
shown to be more generic than its original definition apparently shows. Early work 
showed that the voter model lies at a transition point between a ferromagnetic or- 
dered phase and a paramagnetic disordered phase, such that any small perturbation 
is able to drive it out of the transition and radically change its behavior [26 , 1 81 1231 . 
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Further work on this issue, based on numerical simulations and phenomenologi- 
cal field theories has finally led to the realization that there exists a whole gen- 
eralized voter (GV) universality class encompassing models at an order-disorder 
transition driven by interfacial noise between "dynamically equivalent" absorb- 
ing states. The equivalence can be enforced either by Z2-symmetric local rules, 
or by global conservation of the magnetization II15I I21. The newly observed GV 
universality class is defined in terms of the critical exponents characterizing the 
transition. However, recent research has shown that other characteristics, such as 
the shapes of the exit probability and the consensus times, exhibit nonuniversal 
features, depending on the microscopic details defining the model 11131 . 

The investigation of the existence and the properties of the GV class is based 
on the consideration of nonlinear voter models, in which the flipping probability 
f(x) is a function taking a nonlinear form 01511311 . In the general case, the nonlin- 
earity of the flipping probability can be considered to be modulated by a tunable 
parameter £ and the GV behavior occurs for a critical value £ c , that separates a 
ferromagnetic (ordered) phase, in which the systems orders with a time that scales 
logarithmically with the system size, from a paramagnetic (disordered) phase, in 
which the systems orders at exponentially large times. Few exact results are avail- 
able about these models, since contrary to the standard voter model, nonlinear 
voter models are not in general analytically solvable. Therefore, our understand- 
ing of the GV class is mainly based on numerical simulations and on the analysis 
of field theories, which are studied also numerically or by means of renormal- 
ization group arguments [9|. In this sense, analytical results for nonlinear voter 
models can in general be produced only at the mean field level. 

In this paper we focus in the analysis of the mean field behavior on networks 
of a nonlinear voter model recently introduced in the context of opinion dynamics 
the g-voter model 1121 . This model is characterized by two parameters, represent- 
ing the number of interacting voters q, and a source of noise £. A critical value of 
the parameter £ identifies the GV transition, taking a functional form depending 
on the particular value of q. Considering initially the simplest case of a fully con- 
nected graph, we solve the Fokker-Planck equation that represents the dynamics, 
and obtain analytical information for different values of q. We show that, even 
in the simplest mean-field case, nonlinear voter models can show a very rich be- 
havior, characterized by a complex phase diagram for sufficiently large values of 
q. We then introduce a theoretical treatment, within the heterogeneous mean-field 
framework, for the g-voter dynamics on generic uncorrelated networks. We finally 
turn to the consideration of random regular networks as substrate of the dynamics. 

2 The g-voter model 

On an arbitrary lattice or graph, the g-voter model is defined as follows: Each 
node hosts a spin, with a binary value a, ■ = ± 1 . The system evolves dynamically 
by selecting, at a given time ?, a randomly chosen spin located at, say, node i, with 
state a, . Additionally, q neighbors of node i are also chosen at random, allowing 
for repetition. If all the q neighbors are in the same state, a q , the original spin 
takes the value of the q neighbors, a, := O cl ; otherwise, and with an independent 
probability £.the original spin flips its state, a, := — <7,. In any case, time is updated 
t —¥t + l/N, where N is the total number of sites. 
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It is easy to see that this model is non-linear. Consider the probability / that 
a site flips as a function of the fraction x of neighbors in the opposite state. The 
flipping probability can be easily seen to take the form 

/(jc,e, ? )=x« + e[l-jc*-(l -*)«]. (3) 

For q = 1, we recover, for any value of £, the standard voter model, namely 
f(x, £, 1 ) = x. For q = 2, on the other hand, the flipping probability takes the form 

f(x,e,2) =2ex + (l -2e)x 2 . (4) 

For e = 1 /2, fix, 1/2,2) =x, and therefore linear voter behavior must be observed 
for this combination of parameters in any kind of substrate. In order to observe 
truly nonlinear behavior in the g-voter model (in the sense that f(x) ^ x at the 
transition), we must thus consider values of q larger than 2. 



3 The g-voter model on fully connected graphs 

The case of a fully connected graph corresponds to the classic mean-field behavior, 
in which every node is neighbor of all the other nodes in the system. Therefore, 
choosing q neighbors at random corresponds to the random selection of any q 
nodes in the system. In this case, the state of the system is fully characterized by 
number n of sites in state +1. The behavior of the g-voter model can be derived 
in this regime by means of the exact diffusion formalism developed in Ref. (5). 
Here we follow however a more intuitive approach based on the Fokker-Planck 
equation [19 6 1. The probability of finding the systems in state n at time f fulfills 
the master equation 

= £ T(n\n')P(n',t) - £ T(n'\n)P{n,t), (5) 

where T(n\n') are the transition rates from a state with n' to a state with n positive 
spins. The only non-zero transition rates can be seen to take the form, for a fully 
connected graph, 

T{n+\\n) = j{l-x)f{x,e,q), (6) 

T(n-l\n) = jxf(l-x,e,q), (7) 

T(n\n) = j[\-{\-x)f{x,e,q)-xf{\-x,e,q)}, (8) 

where A = 1 /N is the elementary time scale of the problem and we have defined 
x = n/N. Performing an expansion up to second order in A, the previous master 
equation can be converted into a simplified Fokker-Planck equation, valid in the 
limit of large system sizes, taking the form 
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where we have defined the drift v(x) and the diffusion coefficients D(x), taking 
the forms 

v{x) = (1 -x)f{x,e,q) -xf(l -x,e,q), (10) 
D(x) = i [(1 -x)f(x,e,q) +xf(l -x,e,q)} . (11) 

From the Fokker-Planck equation, application of a standard formalism [19] leads 
to an exit probability E(x) satisfying the differential equation 

v(x)d x E(x) + l -D{x)d 2 x E{x) = 0, (12) 

with boundary conditions E(0) =0 and E(l) = 1, while the average time until 
consensus, T^(x), is given by 

v(x)d x T N (x) + l -D(x)d?T N (x) = -1, (13) 

with boundary conditions T^(0) = 7]y(l) = 0. 

Linear flipping probability f(x) = x implies absence of drift, v(x) =0 and 
then linear exit probability E(x) = x. Moreover, it is easy to see from Eqs. ((TJ) 
and ( [13] ), that a zero drift is necessary and sufficient condition to yield, for any 
diffusion D(x) scaling as 1/N, E(x) =x and 7}y(l/2) ~ N, which are the two 
main signatures of voter behavior at the mean-field level. Let us consider what is 
the behavior of the g-voter model for different values of q. 



3.1 Case q = 2 

For q = 2, the drift and the diffusion coefficient take the form 

v(jc) = -(l-2e)(l-jc)jc(l-2*), (14) 
D(x) = (l + 2e)(l-x)x/N. (15) 

The value £ = 1/2 leads to linear voter behavior, with v(x) = 0. On the other hand, 
the drift takes the form D(x) = 2x(\ —x)/N, leading to a consensus time of the 
form 

T N (x) = -Ar[jdogjt+(l -Jt)log(l -je)] , (16) 

which coincides with the form observed in regular lattices for dimension d > 2, 
Eq. <[TJ. From Eq. ( (14) , we can easily see that £ < 1 /2 implies v(x) < for x < 1/2 
and v(x) > for x > 1/2, which is indicative of a ferromagnetic ordered state. 
The opposite behavior is found for £ > 1/2, indicating a disordered paramagnetic 
phase. Therefore, the point £ = 1/2, corresponding to linear voter behavior, lies at 
a transition between a ferromagnetic and a paramagnetic phase. 



6 




Fig. 1 (color online) Main: Normalized consensus time for the g-voter model on the fully con- 
nected graph for q = 3 at the critical point e = 1/3. For comparison, the entropic form, Eq. (TB) , 
and the one given by Eq. |20} are also plotted. Inset: Normalized consensus time for x = 0.25as 
a function of N, showing the perfect agreement with the prediction of Eq. j20| . Reported results 
from numerical simulations, averaged over 10 4 independent realizations. 



3.2 Case q = 3 



In the case q = 3, we have 

v(x) = (-l+3e)(l-jt)jt(l-2x), 
D(x) = (l-x)x[l+3e-2x(l-x)]/N. 



(17) 
(18) 



The value £ = 1/3 separates again a ferromagnetic phase at £ < 1/3 from a para- 
magnetic phase for £ > 1/3. Exactly at £ = 1/3, we recover v(x) = 0, which leads 
to a linear exit probability, but now D(x) =2(1 -x)x(l -x+x 2 )/N. The diffusion 
term is now different from the case of the linear voter, and as a consequence it 
will lead to a different form of the consensus time as a function of x. The form can 
be recovered by integrating Eq. ( 13 I, upon substituting the diffusion in Eq. (18) , 
namely 



(1 -x)x{\ -x + x z )d{T N (x) = -N. 
The integration of this differential equation takes the form 



(19) 



T N (x) =-N 



xlnx+ (1 — jr)ln(l — x) — 2 m (l ~ x ~ 



1-2* 



tan 



1-2* 



K 



V3 ™" V3 6V3\ (20) 

In Fig.[TJwe plot the normalized consensus time Z#(jc)/2#(1/2) for the q = 3- 
voter model and £ = 1/3, Eq. (20 1, compared with the result from the linear voter 
model, Eq. ( |16| ). The difference induced by the nonlinear form of f(x, 1/3,3) is 
small but noticeable, and indeed it is recovered by means of numerical simulations 
on a fully connected graph (symbols in Fig.[T). 
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3.3 Case q = 4 

In the more interesting case of q = 4, the drift and diffusion coefficient take the 
forms 

v(x) = (l-x>(l-2x)[-l+4e+x(l-x)(l-2e)], (21) 
D{x) = (l-x)jc[l+4e-(3 + 2e>(l-x)]/AT. (22) 

From this form of the drift, it is easy to check that there exists no value of £ that 
allows it to vanish for every value of x. Therefore, there is no region of linear voter 
behavior in the phase diagram of this model. In order to visualize the form of this 
phase diagram, it is convenient to change variables and represent the drift as a 
function of the average magnetization (j) = 2x — 1 . With this transformation, the 
drift takes the form 

vfo) = ~ (1 - </> 2 ) ((1 - 2e)</> 2 + 3 - 14e) . (23) 

A ferromagnetic phase will correspond to the region of values of £ for which 
</>v(</>) > for all </>, while a paramagnetic phase corresponds to <j>v(<j>) < 0. From 
Eq. (j23j, we can see that the nature of the phase is given by the sign of the function 

F e (0) = 3-14£ + (l-2£)</> 2 . (24) 

The ferromagnetic phase corresponds to F e (<j>) > 0, which takes place in the re- 
gion £ < 3/14; the region £ > 1/4 corresponds to a paramagnetic phase with 
F e ((j>) < 0; on the other hand, the region 3/14<£<l/4isa mixed phase, which 
is ferromagnetic or paramagnetic depending on the initial value of </>. In particular, 
initial conditions with magnetization (j) 2 < ^e-l ^ eac ^ t0 a P arama g ne tic behavior, 
while (j) 2 > ^§zrf corresponds to a ferromagnetic phase. 

The ensuing phase diagram for the q = 4-voter model on fully connected 
graphs is depicted, as a function of the density x of + 1 spins, in Fig. [2] Summarizing 
its structure, for £ > 1/4, there is a paramagnetic disordered phase, with T^{x) 
growing exponentially with N and E(x) — > 1/2 in the limit of large N. For £ < 
3/14, on the other hand, a ferromagnetic phase appears, with a fast consensus 
time scaling logarithmically with system size, and an exit probability tending to a 
step-function. Finally, in the region 3/14<£<l/4 there is a mixed phase such 

that for \2x — 1 1 < y ^fp^ i we should expect a paramagnetic disordered phase, 

while for \2x— 1| > y ^jjf > the system should order to consensus given by the 

majority of initial spins. The boundary separating both behaviors in the mixed 
phase is given by an effective e(x) as 



x 2 - x + 1 
2(x 2 -x+2)' 



£(*) = „, 2 , (25) 



such that for £ > e(x), a paramagnetic phase is observed, while £ < e(x) corre- 
sponds to the ferromagnetic phase. 

In Fig.[3]we check the predicted behavior in fully connected networks by mea- 
suring numerically the exit probability E(x) for increasing system size N and three 



Fig. 2 (color online) Phase diagram for the q = 4-voter model on fully connected networks. The 
arrows indicate the direction of the drift for different values of x. 

values of £ representing the three phases discussed above. For £ = 0.30 > 1/4, 
E(x) tends to 1 /2, independent of x, as expected in the paramagnetic phase. For 
£ = 0.15 < 3/14 instead the tendency is toward a step-function &{x— 1 /2), indi- 
cating a ferromagnetic phase. In the intermediate regime (£ = 13/58 = 0.2241 . . .) 
the ferromagnetic phase dominates for x < 1/4 and x > 3/4, while the region 
around x = 1/2 is paramagnetic. 

In analogy with the GV class, one should expect to observe nontrivial behavior 
in the separatrix between the paramagnetic and the ferromagnetic phase. In the 
case of the q = 4-voter model, this separatrix, £ (x) is not a constant, and therefore 
it is impossible to study its behavior in the general case. We shall therefore focus 
on the minimum of the separatrix, given by the points £ = 3/14, and x= 1/2. 
Considering Eq. ( |13| l, and changing again to the total magnetization variable <j> = 
2x— 1, we are led for £ = 3/14 to the equation for the consensus time 



with boundary conditions reading now T((j) = —1) = T((j) = 1) =0, and where 



v(<l>)Tt f (<t>) + ~D(<t>)T%(<j>) 



1 



(26) 



v(</>) 



1 



(27) 



14 
1 



(l-0 2 )(7 + 6</> 2 ) 



(28) 



IN 



By defining 




(29) 



9 





0.5 



0.4 - 



-0.3 



0.2 - 



G-ON = 


20 


B-QN = 


60 


0-ON = 


100 


a-an = 


200 


<~<J N = 


500 


V-VN = 


1000 


t>-t>N = 


2000 




Fig. 3 (color online) Numerical check of the behavior of the q = 4 model on fully connected 
graphs in the different regions of the phase-diagram. The three figures represent the behavior 
of the exit probability for e = 0.30 (top), e = 13/58 (middle) and e = 0.15 (bottom). Different 
curves are for increasing values of the system size N. Results averaged over 10 s independent 
realizations. 
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Fig. 4 (color online) Size dependence of the consensus time for the g-voter model on the com- 
plete graph for x = 1/2, q = 4 and e = 3/14. Numerical simulations of the model (circles) are 
compared with the asymptotic estimate Eq. |3T) (full line). No fitting parameters are introduced. 
Numerical results are averaged over at least 10 4 independent realizations. 



the exact solution to Eq. p6| is given in integral form by 



/(</>) = jdzie-^ £ dz2g{z2)e^\ 



(30) 



We notice that since the first exponential is peaked around zi = as exp ^— , 

it will select values of zi « in the second integral, which can be thus expanded 
for small z\- An estimate of the consensus time can thus be expressed in terms of 
error functions [ 1 1 as follows 



erf 



erf 



M/> 4 
^8" 



exp 



"28 



exp 



N(j) 4 
^8" 



At the point x = 1 /2 (0 = 0), the consensus time takes the form 
T N (x= 1/2) ~ v / 77r7V : erf 



exp 



N 

"28 



(31) 



VritN (32) 



This expression recovers the observed N l l 2 scaling in the large limit [ 12 1, and 
is found in good agreement with numerical simulation results, as shown in Fig. [4] 
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Fig. 5 (color online) Left: Consensus time for the q- voter model (q = 4, e = 3/14) on the 
complete graph as a function of the initial condition x for different system sizes. The results of 
a numerical integration of the exact formula Eq. < f30j > are plotted for increasing values of the 
system size. Right. Consensus time obtained from numerical simulations as a function of N 
for the initial condition x = 1/3. Simulation results are averaged over at least 10 4 independent 
realizations. The hump in the Tn function is indicative a crossover effect. 



For x 7^ 1 /2 the system is in the ferromagnetic phase and one expects the con- 
sensus time to scale logarithmically for large values of the system size N. By direct 
numerical integration of Eq. ( |30| l, the value of the consensus time for any value 
of x and N can be computed, see Fig. [5] (left panel). The crossover between the 
behavior at x = 1/2, 7Jv ~ \/N, and the pure logarithmic growth which asymptoti- 
cally prevails for all x 7^ 1 /2, give rise to a nonmonotonicity in the functional form 
of 7/v as a function of N, see Fig. [5] (right panel). This is a curious phenomenon: 
there is a short range of values of N for which the consensus time decreases as N 
is increased. 

Concerning the exit probability, the solution of the corresponding equation 
leads to the form, again expressed in terms of the magnetization <j>, 

= „o.i _n 2 ,_ — : . < 33 ) 



ji;e-T2^(7 + 6z 2 ) 1N / 12 d. 



z 



Expanding the integral in powers of 1 /N, and keeping only the first order, we are 
led to 

J+l e-^l^dz 



£ W - r +i , ■ (34) 



Focusing on the case < (x < 1/2), due to the symmetry of the expression, the 
integral can be performed, yielding 



I 1 U' 28 J 1 V4'28y 



E ($) - o wn „/i ^ (35) 



2 r (i)- r U>2S, 
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Fig. 6 (color online) Exit probability for the q- voter model (for q = 4) at the transition point 
£ = 3/14 on fully connected networks. Lines report the theoretical prediction Eq. (35) . Symbols 
are the results of numerical simulations averaged over 10 4 realizations of the dynamics. 



where F (z) is the Gamma function and F (a ,z) is the incomplete Gamma function 
|Q~). In Fig. [6] we plot Eq. ( |35[ ) together with numerical simulations, displaying 
again a good agreement for values of N larger than 100. 



3.4 Phase diagram for general q 

Following the same arguments as in the previous section, we can sketch a phase 
diagram for the general case of the g-voter model by considering the sign of the 
function <j>v(<j>). Considering the form of the drift in Eq. (jTOj and the flipping 
probability in Eq. simple calculations allow to show that the quantity <j>v(<j>) 
is positive for £ < £, with £ taking the form (expressed as a function of x for 
simplicity) 



e(x) 



;i-2*)[l -*<?-(! -*)<?]' 



(36) 



In Fig.|7]we plot the shape of the phase boundary e[x), Eq. ( |36| ), as a function of 
x, for different values of q. We observe that 



lim e(x) 



1 

J 

9 



lim e(x) 

*->l/2 



Or 
21- 



(37) 



We thus obtain that, for the general q case, the g-voter model is in a ferromag- 
netic phase for e < , and in a paramagnetic phase for e > | . In the region 

J^y < £ < |, the system is in a mixed phase whose nature depends on the ini- 
tial conditions. In the limit of large q, the boundary of region for ferromagnetic 
behavior decreases exponentially, while the boundary for paramagnetic behavior 
decreases only algebraically. 
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Fig. 7 (color online) Plot of the phase boundary e(x), Eq. {36}, as a function of x, for different 
values of q. 



4 The q- voter model on networks: General heterogeneous mean-field theory 

The analytical treatment of voter-like models on complex networks, beyond the 
simplest fully connected graph, is traditionally based on the heterogeneous mean- 
field (HMF) approach, a powerful tool for the theoretical analysis of general dy- 
namical processes on heterogeneous substrates [4 16|. HMF theory is based on a 
fundamental approximation: The real (quenched) network is coarse-grained into 
an annealed one [16|, which disregards the specific connection pattern and postu- 
lates that the class of degree k is connected to the class of degree k! with condi- 
tional probability P(k'\k) Q. 

The application of the HMF formalism to the g-voter model follows from an 
extension of the case of standard voter model II291I28II . taking into account the 
intrinsic multipoint nature of the g-voter dynamical rules. Indeed, the annealed 
network approximation is based on the assumption that the network gets com- 
pletely rewired between any two dynamical time steps. In dynamics in which a 
vertex interacts with a single neighbor, or with all of them, the numerical (and 
analytical) implementation goes by selecting one or k, (according to the case) 
neighbors, randomly chosen with probability kjP(kj)/ (k), where we are for sim- 
plicity making the further assumption that the network lacks degree correlations, 



In the g-voter model, we choose q neighbors, possibly with repetition, among 
the k{ neighbors of a vertex. If we choose these q neighbors at random among all 
the vertices in the network, we are implicitly assuming that there is rewiring of the 



i.e. P(k'\k)=k'P(k')/(k) El- 
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network even during the duration of an elementary dynamic step and this under- 
estimates the probability that repetitions occur, i.e. a neighbor is selected multiple 
times. The proper way to eliminate this inconsistency is to select at random a fixed 
set of ki neighbors among all the vertices in the network, and then to pick up the q 
participants to the dynamic step only from the fixed set of kj neighbors previously 
selected. 

Let us thus consider a generic sparse network. Assuming no correlations and 
annealed approximation, the probability that a randomly chosen neighbor of a 
vertex is in state +1 is given by 

e = E^ft, (38) 

where x% is the probability that a vertex of degree k is in state +1. In the set 
of k neighbors of a vertex of degree k, chosen to make an elementary update of 
this vertex, the probability that we choose n in state +1 is given by the binomial 
distribution B kM (6) = (*) 0"(1 - 9) k ~". In this case, if vertex k is in state -1, it 
will flip with probability 

fk(n,q) = (jj +£ 

while if it is in state +l, it will flip with probability f k (k — n,q). Therefore, the 
probabilities TI(k, a) that a microscopic update will flip a spin of degree k and 
state s [24 1 are given by the elementary flipping probabilities, fk(n,q) and fk(k — 
n,q), averaged over the binomial distribution B k „(6), and take the form 

n(k,+i)=p(k)x k J £ ( k )e n (i-e) k -"f k (n, q ), (40) 

n=Q W 

n(k-l)=P(k)(l-x k )Y,( k )d n (l-e) k - n f k (k-n,q), (41) 

where the first terms stand for the probability of selecting a vertex to update of 
degree k, P(k), and the probability that it is in state +1 (jc*), or —1 (1 — x^), re- 
spectively. The HMF equation for the time evolution of the density x k takes the 
form 

x k = ^-[n(k,-i)-n(k,+i)}. (42) 

From this solution, one can find information about the exit probability by look- 
ing for conservation laws of the dynamics in Eq. ( [42) 1, while the consensus time 
follows from an immediate one-step calculation complemented with an adiabatic 
approximation, which allows to replace the multiple degrees of freedom x k by 
single slowly varying conserved quantity ['29,28). 

For q = 2 and any network, the application of this formalism leads to the trivial 
result presented in Sec. [2] For values of q > 2, the application of this formalism 
is hampered by the fact that the (possible) conserved quantities in Eq. ( |4"2"| ) are not 
linear functions of x k , which implies that the transformations leading to the equa- 
tion for the consensus time cannot be applied directly [24|. We therefore consider 
in the following section the simpler case of random regular networks, where a full 
analytical solution can be obtained for q = 3. 



/n\i (k — n\ 
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5 The g-voter model on random regular networks 

Random regular networks (RRN) are networks where all nodes have exactly the 
same degree k, while links are randomly distributed among them, avoiding self- 
connections and multiple connections. In this case, we have P(k) = 8i<.k, where K 
is the average (constant) degree of the network. In this particular network, since 
all vertices share the same value of the degree, a great simplification arises, as we 
can set xk = 6 = x, the total average density of vertices in state +1 . Additionally, 
since the density is homogeneous, we can sidestep the general formalism from 
HMF theory and apply directly the Fokker-Planck approach presented in Sec. [3] 

In the case q = 2, the general considerations from Sec.|4]indicate the presence 
of linear voter behavior for £ = 1/2. 



5.1 Case q = 3 



In the more interesting case of q = 3, we are led, from Eq. 
a drift and a diffusion coefficient of the form 



v(x) = 
D(x) = 



([TO]), and ([TT), to 

(43) 



;i -x)x(l-2x){K-l)[(3e-l)K-\] 
K 2 ' 
(1 -x)x [3e(K-l)K + K 2 -2(K-2)(K- l)x(l -x) + l] 

Yn 



The drift turns out now to vanish for a value £ c that depends on K as 

, * K+l 
e c (K) 



3K 



(44) 



(45) 



This value tends to 1 /3, corresponding to a fully connected network, in the limit 
K — > oo, but is always larger than 1/3 for any finite K. This observation indicates 
that the structure of a RRN induces more order than a fully connected topology, 
and thus that a larger disordering parameter £ is needed to cross over to the disor- 
dered paramagnetic phase. 

For this critical value the diffusion coefficient is 



D(x) 



2(1- x)x [K 2 - {K - 2) (K - 1 )x{ 1 - x)] 
" K 2 N " 



Using Eq. ( |13| l we are led to the equation for the consensus time 

(1 -x)x[\ - ax(l -x)) d*T N (x) = —N, 
where a = (K — 2) (K — l)/K 2 , whose solution is 



(46) 



T N (x) = -N 



jclnx+ (1 — x)ln(l —x) — -ln(l — ax+x 2 a) 



(l-2x)^tan-'(( 1 =gy) 



a tan 



l-a J 



(47) 
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Fig. 8 (color online) Exit probability E(x = 0.25) for the g-voter model with q = 3 on a random 
regular network with K = 4 as a function of N for different values of e. The results are obtained 
by averaging over 10 5 numerical simulations of the dynamics. The plateau singles out the critical 
value e c . 



Again, we obtain a non-entropic form, modulated by the factor a, and that yields 
the mean-field result Eq. ( pO] l in the limit a — )■ 1 (K — >• oo). 

To check the above results we have simulated the q = 3 -voter model on RRN 
with K = 4. In this case, HMF theory predicts a transition at 



12 



0.4166666. 



(48) 



Numerical results lead however at a different, much lower, value. In order 
to estimate numerically the value of E c we have followed the following approach, 
based on the behavior of the exit probability E (x) 113L Indeed, £ > £ c corresponds 
to a disordered paramagnetic phase with, for asymptotically large systems, E(x) = 
1/2, while £ < £ c corresponds to an ordered ferromagnetic phase, where E(x) = 
®(x), the Heaviside theta function. Therefore, focusing on an initial density x < 
1 /2, we should observe E(x) — > 1/2 for £ > e c , E(x) — > for £ < £ c , and E(x) — > 
const < 1/2 for £ = £ c when increasing the system size L. In Fig. [8] we report 
the exit probability for x = 0.25 and different values of £ as a function of the 
lattice size ,/V. A plateau is obtained for £ ~ 0.3495, while larger (smaller) values 
of £ lead to an increase (decrease) of E(x = 0.25) with ,/V. We conclude that the 
critical point is in this case located at £ c = 0.3495 ±0.0005, where the error bars 
are obtained from the values of £ leading to an increasing or decreasing behavior 
for E(x = 0.25). Fig.|8]also indicates that E(x = 0.25) < 0.25, i.e. the global form 
of the exit probability at the transition is not linear, in agreement with the findings 
in d — 2 1131 . Not unexpectedly, also the dependence of the consensus time on x 
(Fig. [9) turns out to be different from the analytical prediction, Eq. (47 1). 
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Fig. 9 (color online) Main: Normalized consensus time for the q- voter model with q = 3 on a 
random regular network with K = 4, e = e c = 0.3495 and different values of N. The numerical 
results are averaged over 10 5 realizations of the dynamics. Inset: Normalized consensus time for 
x = 0.25 as a function of N, showing the disagreement with the theoretical predictions. 



5.2 Case q = A 

The case of q = 4 in RRN leads to rather complex expressions. In particular, as a 
function of the magnetization <j>, the drift takes the form 

(K-l)(j> (</> 2 -l) (UK 2 e~3K 2 + (K-3)(K-2)(2e-l)(j) 2 + 2Ke-9K-4e + 2) 

VW = ' 16F 

(49) 

As in the case of the fully-connected graph, no value of £ can make this expression 
zero for any <j>. Therefore we exclude again for this value of q the presence of a 



linear voter point. Performing a similar analysis as in Sec. 3.3 we can however 
determine the phase diagram by looking at the sign of the function </>v(</>). It turns 
out that </>v(</>) is positive for £ smaller than the threshold £, which, as a function 
of x, takes the form 

1 {K-2)K 
eW_ 2~ 2{K 2 ({x-\)x + 2)+K{-5{x-\)x-\)+6{x-\)x+\y { ' 

The phase-diagram is perfectly analogous to the one illustrated in Fig. |2j with a 
ferromagnetic phase for 

,. , 3K 2 +9K-2 
£< <% £{x) = 2(lK 2 + K-2y (51) 

a paramagnetic phase for 



(52) 
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Fig. 10 (color online) Exit probability of the g-voter model for q = 4 on random regular net- 
works with K = 4, slightly above (e = 0.25,left) and below (£ = 0.24, right) the transition. 
Different symbols are for increasing values of the system size N. It turns out that the system is 
in the ferromagnetic phase for any x for e = 0.24 and in the paramagnetic phase for any x for 
e = 0.25. The numerical results are averages over 10 5 realizations of the dynamics. 



and an intermediate mixed phase, whose nature depends on the initial condition x. 

The interesting question is whether the real g-voter dynamics in quenched 
RRN networks obeys (at least qualitatively) these theoretical predictions, in par- 
ticular with respect to the existence of the mixed phase. We have checked this by 
means of numerical simulations in RRN with K = 4. In this case, Eq. ( 50 1 predicts 
a mixed phase for values of £ in the range 21/58 < £ < 41 /l 14. Numerically we 
observe instead (Fig. 10 1 that there is a simple transition (for £ c « 0.247) separat- 
ing a purely ferromagnetic phase for £ < e c (where E(x) tends to a step-function 
@{x— 1 /2)), from a purely paramagnetic phase for £ > £ c (where E(x) tends to 
1/2 for any x). Hence we conclude that the fluctuations introduced by the RRN 
quenched topology are sufficient to make the behavior of the g-voter model similar 
to what happens on finite-dimensional lattices [1211131 and qualitatively different 
from what happens on fully-connected graphs. We expect similar conclusions to 
hold for any q > 4 and for any K < N — 1 . 



6 Conclusions 

In this paper we have presented a detailed mean-field analysis of a model for 
social dynamics belonging to the generalized voter universality class. A rich phe- 
nomenology emerges, in particular for q > 4, with a mixed phase separating para- 
magnetic and ferromagnetic regions of the phase diagram. When the dynamics 
occurs on networks a natural question arising is whether the behavior is of mean- 
field finite dimensional type. From the analysis of random regular networks it 
turns out that the minimal disorder present in those graphs is sufficient to de- 
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stroy the mixed phase for q > 4 and restore a behavior qualitatively similar to 
what happens in lattices, with a single transition point separating ordered and dis- 
ordered regions of the phase diagram. Also in the case q = 3, the dynamics on 
RRN is more akin to finite dimensional lattices (nonlinear exit probability at the 
transition) than complete graphs (where E(x) = x). Questions about the behavior 
on complex, heterogeneous networks remain open. In that case it is possible that 
more interesting phenomena could arise, due to the interplay between the value of 
q and the degree k, which has large fluctuations from node to node. 
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